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We report on flux confinement effects in superconducting submicron line, loop and 
dot structures. The main idea of our study was to vary the boundary conditions 
for confinement of the superconducting condensate by taking samples of different 
topology and, through that, modifying the lowest Landau level Elll(H). Since 
the critical temperature versus applied magnetic field T C {H) is, in fact, Elll(H) 
measured in temperature units, it is varied as well when the sample topology is 
changed. We demonstrate that in all studied submicron structures the shape of the 
T C (H) phase boundary is determined by the confinement topology in a unique way. 
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1. Introduction 

Recent impressive progress in nanostructuring (e-beam lithography, single atom manipulation 
with a scanning tunneling microscope tip, etc.) has made it possible to control the quantization 
effects in nanostructures by varying their size and topology. The most commonly used so far for 
this purpose are semiconducting structures, where an elegant approach developed by Landauer JO 
relates directly mesoscopic transport to the quantum transition probability, which is determined 
by the specific configuration of quantum levels and available tunneling barriers. Changing the con- 
finement potential via nanostructuring allows to tune both the quantum levels and the tunneling 
probabilities. Extended by Biittiker ^| to the multi-terminal measurement geometry, the Landauer- 
Buttiker formalism has been widely and very successfully used in the interpretation of numerous 
mesoscopic experiments 

In comparison to the semiconducting and normal systems, superconducting structures have been 
studied much less so far. In this report, it is worth first to ask a few simple questions like: why do 
we want to make such structures, what interesting new physics do we expect, and why do we want 
to focus on superconducting (and not, for example, normal metallic) nanostructured materials? 

First of all, by making low dimensional systems, one creates an artificial potential in which charge 
carriers or flux lines are confined. The confinement length scale La of an elementary "plaquette" 
A, gives roughly the expected energy scale E = h 2 Tr 2 n 2 /(2mL A ). The concentration of charge 
carriers or flux lines can be controlled by varying the gate voltage in two-dimensional (2D) electron 
gas systems [Q or the applied magnetic field in superconductors [|). In this situation, different 
commensurability effects between the fixed number of elements A in an array and a tunable number 
of charge or flux carriers are observed. 

Secondly, modifying the sample topology in those systems creates a unique possibility to impose 
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the desired boundary conditions, and thus almost "impose" the properties of the sample. A Fermi 
liquid or a superconducting condensate confined within such a system will be subjected to severe 
constraints and, as a result, the properties of these systems will be strongly affected by the boundary 
conditions. 

While a normal metallic system should be considered quantum-mechanically by solving the 
Schrodinger equation: 

-ih\7 - eA) 2 ^ + U ^ = E ^ , (f) 



J-f- 
2m V 



a superconducting system is described by the two coupled Ginzburg-Landau (GL) equations: 

1 -(-zfiV -e*I) 2 * s + /3|# s | 2 * s = -a$> s (2) 



2m 



j s = V x h = 



2m* 



V*(-ihV - eM)* s + $ s (iW - e*i)f* , (3) 



with A the vector potential which corresponds to the microscopic field h — rotA/no, U the potential 
energy, E the total energy, a a temperature dependent parameter changing sign from a > to a < 
as T is decreased through T c , (3 a positive temperature independent constant, m* the effective mass 
which can be chosen arbitrarily and is generally taken as twice the free electron mass m. 

Note that the first GL equation (Eq. (0)), with the nonlinear term /3|\E r s | 2, J' s neglected, is the 
analogue of the Schrodinger equation (Eq. TO) with {7 = 0, when making the following substitutions: 
4 r s (4 r ), e*(e), —a{E) and m*(m). The superconducting order parameter SSf s corresponds to the wave 
function the effective charge e* in the GL equations is 2e, i.e. the charge of a Cooper pair; the 
temperature dependent GL parameter a 

-a= (4) 

plays the role of E in the Schrodinger equation. Here £(T) is the temperature dependent coherence 
length: 

T \ -1/2 



£(T)=i(0) j . (5) 

The boundary conditions for interfaces between normal metal-vacuum and superconductor- 
vacuum are, however, different (Fig. Q): 

= (6) 
{-ihV - e*A)^ s = (7) 

_L,6 

i.e. for normal metallic systems the density is zero, while for superconducting systems, the gradient 
of (for the case ^4 = 0) has no component perpendicular to the boundary. As a consequence, 
the supercurrent cannot flow through the boundary. The nucleation of the superconducting con- 
densate is favored at the superconductor/ vacuum interfaces, thus leading to the appearance of 
superconductivity in a surface sheet with a thickness £(T) at the third critical field H c3 (T). 

For bulk superconductors the surface-to- volume ratio is negligible and therefore superconductivity 
in the bulk is not affected by a thin superconducting surface layer. For submicron superconductors 
with antidot arrays, however, the boundary conditions (Eq. (^)) and the surface superconductivity 
introduced through them, become very important if La < £,{T). The advantage of superconducting 
materials in this case is that it is not even necessary to go down to the nanometer scale (like for 
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Fig. 1. Boundary conditions for interfaces between normal metal-vacuum and superconductor-vacuum. 



normal metals), since for La of the order of 0.1-1.0 fim the temperature range where La < £,(T), 
spreads over 0.01 - 0.1 K below T c due to the divergence of £(T) at T -> T c0 (Eq. (§)). 

In principle, the mesoscopic regime La < C(^) can De reached even in bulk superconducting 
samples with La ~ 1 cm, since £(T) diverges. However, the temperature window where La < £(^) 
is so narrow, not more than ~ 1 nK below T cu , that one needs ideal sample homogeneity and a 
perfect temperature stability. 

In the mesoscopic regime, La < £{T), which is easily realized in (perforated) nanostructured ma- 
terials, the surface superconductivity can cover the whole available space occupied by the material, 
thus spreading superconductivity all over the sample. It is then evident that in this case surface 
effects play the role of bulk effects. 

Using the similarity between the linearized GL equation (Eq. (g)) and the Schrodinger equation 
(Eq. ([!])), we can formalize our approach as follows: since the parameter -a (Eqs. (||) and (^)) plays 
the role of energy E (Eq. ([l])), then the highest possible temperature T C (H) for the nucleation of 
the superconducting state in the presence of a magnetic field H always corresponds to the lowest 
Landau level Elll(H) found by solving the Schrodinger equation (Eq. (Q)) with "superconducting" 
boundary conditions (Eq. (0)). 

Figure || illustrates the application of this rule to the calculation of the upper critical field H C 2 (T) : 
indeed, if we take the well-known classical Landau solution for the lowest level in bulk samples 
Elll(H) — hu/2, where u) = e* n^H/m* is the cyclotron frequency. Then, from -a — Elll(H) we 
have 



and with the help of Eq. 



2m* £ 2 (T) 
we obtain 

fJ,oH c2 (T) 



Tilo 
~~2 



H=H c2 



4> 



2< 2 (T) 



(8) 



(9) 



with <I>o = h/e* — h/2e the superconducting flux quantum. 

In nanostructured superconductors, where the boundary conditions (Eq. (|7|)) strongly influence 
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Fig. 2. Landau level scheme for a particle in a magnetic field. From the lowest Landau level Elll(H) the second 
critical field H C 2(T) is derived (solid line). 

the Landau level scheme, Elll(H) has to be calculated for each different confinement geometry. 
By measuring the shift of the critical temperature T C (H) in a magnetic field, we can compare the 
experimental T C (H) with the calculated level Elll(H) and thus check the effect of the confinement 
topology on the superconducting phase boundary for a series of nanostructured superconducting 
samples. The transition between normal and superconducting states is usually very sharp and there- 
fore the lowest Landau level can be easily traced as a function of applied magnetic field. Except when 
stated explicitly, we have taken the midpoint of the resistive transition from the superconducting 
to the normal state, as the criterion to determine T C (H). 

2. Flux confinement in individual structures: line, loop and dot 

In this section we present the experimental T C {H) phase boundary measured in superconducting 
aluminum mesoscopic structures with different topologies with the same width of the lines (w — 
0.15 fim) and film thickness (t = 25 nm). The magnetic field H is always applied perpendicular to 
the structures. 

2.1. Line structure 

In Fig. ||a the phase boundary T C (H) of a mesoscopic line is shown. The solid line gives the T C (H) 
calculated from the well-known formula [pi : 




(10) 



which, in fact, describes the parabolic shape of T C (H) for a thin film of thickness w in a parallel 
magnetic field. Since the cross-section, exposed to the applied magnetic field, is the same for a film 
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of thickness w in a parallel magnetic field and for a mesoscopic line of width w in a perpendicular 
field, the same formula can be used for both configurations Q. Indeed, the solid line in Fig ||a is a 
parabolic fit of the experimental data with Eq. ( |To| ) where £(0) = 110 nm was obtained as a fitting 
parameter. The coherence length obtained using this method, coincides reasonably well with the 
dirty limit value £(0) = 0.85(£ ^) 1/2 = 132 nm calculated from the known BCS coherence length 
£o = 1600 nm for bulk Al j^] and the mean free path I = 15 nm, estimated from the normal state 
resistivity p at 4.2 K [§. 

We can use also another simple argument to explain the parabolic relation T C (H) oc H 2 , since the 



expansion of the energy E(H) in powers of H, as given by the perturbation theory, is 10|: 

E(H) = E Q +A 1 LH + A 2 SH 2 + ■■■ (11) 

where A\ and A2 are constant coefficients, the first term Eq represents the energy levels in zero 
field, the second term is the linear field splitting with the orbital quantum number L and the third 
term is the diamagnetic shift with S, being the area exposed to the applied field. 




Fig. 3. The measured superconducting/normal phase boundary as a function of the reduced temperature T C LH) /T c o 
for a) the line structure, and b) the loop and dot structure. The solid line in (a) is calculated using Eq. (Iicj) with 
£(0) = 110 nm as a fitting parameter. The dashed line represents T C (H) for bulk Al. 



For the topology of the line with a width w much smaller than the Larmor radius rjj 3> w, any 
orbital motion is impossible due to the constraints imposed by the boundaries onto the electrons 
inside the line. Therefore, in this particular case L = and E(H) = E + A 2 SH 2 , which immediately 
leads to the parabolic relation T c cx H 2 . This diamagnetic shift of T C (H) can be understood in terms 
of a partial screening of the magnetic field H due to the non-zero width of the line [jllj . 

2.2. Loop structure 

The T C {H) of the mesoscopic loop, shown in Fig. |^b, demonstrates very distinct Little-Parks (LP) 
oscillations Jl2[ superimposed on a monotonic background. A closer investigation leads to the con- 
clusion that this background is very well described by the same parabolic dependence as the one 
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which we just discussed for the mesoscopic line [Q (see the solid line in Fig. ||a). As long as the width 
of the strips w, forming the loop, is much smaller than the loop size, the total shift of T C (H) can be 
written as the sum of an oscillatory part and the monotonic background given by Eq. ( [To| ) ^ |l3| : 



T C (H) = T d 



1 



3 V $o / R 2 V $ o 



(12) 



where R 2 = R\ R2 is the product of inner and outer loop radius, and the magnetic flux threading 
the loop $ = ttR 2 ijlqH . The integer n has to be chosen so as to maximize T C (H) or, in other words, 
selecting E LLL {H). 

The LP oscillations originate from the fluxoid quantization requirement, which states that the 
complex order parameter ^ s should be a single-valued function when integrating along a closed 
contour 

Vcp-dl = n2n n = • • • , -2, -1, 0, 1, 2, ■ • ■ (13) 

where we have introduced the order parameter \& s = \& B \ exp {up). Fluxoid quantization gives rise 
to a circulating supercurrent in the loop when $ ^ n<S>Q, which is periodic with the applied flux 
$/$ . 

Using the sample dimensions and the value for £(0) obtained for the mesoscopic line (with the 
same width w = 0.15 /J,m), the T C (H) for the loop can be calculated from Eq. ( |l2| ) without any free 
parameter. As shown in Fig. |^b, the agreement with the experimental data is very good. 

Another interesting feature of the mesoscopic loop or other structures is the unique possibility 
they offer for studying nonlocal effects JTij. In fact, a single loop can be considered as a 2D artificial 
quantum orbit with a fixed radius, in contrast to Bohr's description of atomic orbitals. In the 
latter case the stable radii are found from the quasiclassical quantization rule, stating that only 
an integer number of wavelengths can be set along the circumference of the allowed orbits. For a 
superconducting loop, however, supercurrents must flow, in order to fulfill the fluxoid quantization 
requirement (Eq. (|13|)), thus causing oscillations of T c versus H. 

In order to measure the resistance of a mesoscopic loop, electrical contacts have, of course, to 
be attached to it, and as a consequence the confinement geometry is changed. This "disturbing" or 
"invasive" aspect can now be exploited for the study of nonlocal effects [B. Due to the divergence 
of the coherence length f(T) at T = T c0 (Eq. (|)) the coupling of the loop with the attached leads 
is expected to be very strong for T — > T c q. 

Fig. |J shows the results of these measurements. Both "local" (potential probes V1/V2 across the 
loop) and "nonlocal" (potential probes V1/V3 or V2/V4 aside of the loop) LP oscillations are clearly 
observed. For the "local" probes there is an unexpected and pronounced increase of the oscillation 
amplitude with increasing field, in disagreement with previous measurements on Al microcylin- 
ders |l^| . In contrast to this, for the "nonlocal" LP effect the oscillations rapidly vanish when the 
magnetic field is increased. 

When increasing the field, the background suppression of T c (Eq. @) results in a decrease of 
£(T). Hence, the change of the oscillation amplitude with H is directly related to the temperature- 
dependent coherence length. As long as the coherence of the superconducting condensate extends 
over the nonlocal voltage probes, the nonlocal LP oscillations can be observed. 

The importance of an " arm" attached to a mesoscopic loop was already demonstrated theoretically 
by de Gennes in 1981 |lf|. For a perfect ID loop (vanishing width of the strips) adding an "arm" 
will result in a decrease of the LP oscillation amplitude, what we indeed observed at low magnetic 
fields, where £(T) is still large. With these experiments, we have proved that adding probes to a 
structure considerably changes both the confinement topology and the phase boundary T C (H). 
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T o(K) 

Fig. 4. Local (V1/V2) and nonlocal phase boundary (V1/V3 or V2/V4) measurements. The transport current flows 
through The solid and dashed lines correspond to the theoretical T C (H) of an isolated loop and a one- 

dimensional line, respectively. The inset shows a schematic of the mesoscopic loop with various contacts (P = OAfim). 



The effect of topology on T C (H), related to the presence of the sharp corners in a square loop, 
has been considered by Fomin et al. [fl6| , |l7|| . In the vicinity of the corners the superconducting 
condensate sustains a higher applied magnetic field, since at these locations the superfluid velocity is 
reduced, in comparison with the ring. Consequently, in a field-cooled experiment, superconductivity 
will nucleate first around the corners []l6|| . Eventually, for a square loop, the introduction of a local 
superconducting transition temperature seems to be needed. As a result of the presence of the 
corner, the H c s{T) of a wedge with an angle 9 [H will be strongly enhanced at the corner resulting 
in the ratio H c3 /H c2 « 3.79 for 9 « 0.44 tt §J]. 

2.3. Dot structure 

The Landau level scheme for a cylindrical dot with "superconducting" boundary conditions 
(Eq. (0)) is presented in Fig. |^. Each level is characterized by a certain orbital quantum num- 
ber L where ^ s = \^ s \ exp (=fiL<p) |l9| . The levels, corresponding to the sign "+" in the argument 
of the exponent are not shown since they are situated at energies higher than the ones with the 
sign "-" . The lowest Landau level in Fig. |H| represents a cusp-like envelope, switching between differ- 
ent L values with changing magnetic field. Following our main guideline that Elll(H) determines 
T C {H), we expect for the dot the cusp- like superconducting phase boundary with nearly perfect 
linear background. The measured phase boundary T C (H), shown in Fig. |^b, can be nicely fitted 
by the calculated one (Fig. 0)^ thus proving that T C (H) of a superconducting dot indeed consists 
of cusps with different L'b pC|. Each fixed L describes a giant vortex state which carries L flux 
quanta <5>q. The linear background of the T C (H) dependence is very close to the third critical field 



H c z(T) ~ 1.69 H C 2{T) |21 . Contrary to the loop, where the LP oscillations are perfectly peri- 
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Fig. 5. Energy level scheme versus normalized flux <&/&q for a superconducting cylinder in a magnetic field parallel 
to the axis. The cusp-like H C ^(T) line is formed due to the change of the orbital quantum number L. 



The lower critical field of a cylindrical dot H^° l corresponds to the change of the orbital quantum 
number from L = to L = 1, i.e. to the penetration of the first flux line [p3[ : 

lioHg* = 1.924-^. (14) 

For a long mesoscopic cylinder described above, demagnetization effects can be neglected. On 
the contrary, for a thin superconducting disk, these effects are quite essential |^4[ ^5], ^6|. For a 
mesoscopic disk, made of a Type-I superconductor, the phase transition between the superconducting 
and the normal state is of the second order if the expulsion of the magnetic field from the disk can 
be neglected, i.e. when the disk thickness is comparable with £ and A. When the disk thickness 
is larger than a certain critical value first order phase transitions should occur. The latter has 
been confirmed in ballistic Hall magnetometry experiments on individual Al disks [^7], [2l| , A 
series of first order transitions between states with different orbital quantum numbers L have been 



seen in magnetization curves M(H) 27 in the field range corresponding to the crossover between 
the Meissner and the normal states. Besides the cusplike H c3 (T) line, found earlier in transport 
measurements p0| , transitions between the L — 2 and L = 1 states have been observed 27 by 



probing the superconducting state below the T C (H) line with Hall micromagnetometry. Still deeper 
in the superconducting area the recovery of the normal <!>o-vortices and the decay of the giant vortex 
state might be expected (^6|. The former has been considered in Ref. in the London limit, by 
using the image method. Magnetization and stable vortex configurations have been recently analyzed 



in mesoscopic disks in Refs. 24, E5l Uu 
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3. Conclusions 

We have carried out a systematic study of and quantization phenomena in submicron structures of 
superconductors. The main idea of this study was to vary the boundary conditions for confining the 
superconducting condensate by taking samples of different topology and, through that, to modify 
the lowest Landau level E LLL (H) and therefore the critical temperature T C {H). Different types 
of individual nanostructures were used: line, loop and dot structures. We have shown that in all 
these structures, the phase boundary T C (H) changes dramatically when the confinement topology 
for the superconducting condensate is varied. The induced T C (H) variation is very well described by 
the calculations of Elll{H) taking into account the imposed boundary conditions. These results 
convincingly demonstrate that the phase boundary in T C {H) of mesoscopic superconductors dif- 
fers drastically from that of corresponding bulk materials. Moreover, since, for a known geometry 
Elll{H) can be calculated a priori, the superconducting critical parameters, i.e. T C (H), can be con- 
trolled by designing a proper confinement geometry. While the optimization of the superconducting 
critical parameters has been done mostly by looking for different materials, we now have a unique 
alternative - to improve the superconducting critical parameters of the same material through the 
optimization of the confinement topology for the superconducting condensate and for the penetrating 
magnetic flux. 
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